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Electromagnetic  Instability  of  an  Intense  Beam 
in  a  Quadrupole  Focusing  System 


Cha-Mei  Tang,  Thomas  Suiyden’and  Jonathan  Krall 
Plasma  Physics  Division,  Naval  Research  Laboratory 
Washington,  DC  20375-5000 


Abstract 

Discrete  quadrupole  focusing  systems  are  subject  to  the 
electromagnetic  instability  previously  found  in  helical 
quadrupole  strong  focusing  systems.  The  dispersion  re¬ 
lation  is  derived.  Analytical  growth  rates  are  obtained. 
Instability  boundaries  are  established.  Numerical  solu¬ 
tions  of  the  dispersion  relation  results  show  that  the  overall 
stability  properties  are  not  favorable  for  long  pulse,  high- 
current  electron  beams. 

I.  Dispersion  Relation 

Discrete  quadrupole  focusing  systems,  also  called 
FODO  lattices,  have  been  used  to  transport  charged  parti¬ 
cle  beams  for  a  variety  of  applications.  Helical  quadrupole 
(stellarator)  focusing  systems  are  subject  to  an  electromag¬ 
netic  instability,  which  we  referred  to  as  the  three-wave 
instability.1-3  This  has  been  observed  experimentally.4 
Here,  we  consider  the  corresponding  instability5  for  trans¬ 
verse  perturbations  of  a  beam  centroid  interacting  with  a 
FODO  lattice  field  and  a  TE\\  waveguide  mode. 

We  consider  an  alternating  gradient  quadrupole  field 
(BqT,Bqy),  where 

Bqz  —  — Bqkqf(z)y ,  Bqy  =  —Bqkqf(z)x,  (la  —  b) 

Bqkq  is  the  peak  quadrupole  field,  f(z )  is  periodic  with 
maximum  value  of  one,  kq  =  2ir/Xq  and  A,  is  the  pe¬ 
riod  of  the  quadrupole  field.  The  representation  for  the 
quadrupole  field  in  Eqs.  (la-b)  is  valid  near  the  z-axis, 
i.e. ,  (z2  +  y2)1/2  <<  A9/2tt.  In  equilibrium,  the  electron 
beam  travels  along  the  axis  of  a  circular  waveguide  at  ve¬ 
locity  va  and  is  monoenergetic  with  7„  =  (1  —  /32)-1/2, 
where  0„  —  v„/c.  For  simplicity,  we  let  /(z)  =  cos (nkqz). 

W<*  assume  the  electron  beam  propagates  within  a  per¬ 
fectly  conducting  cylindrical  waveguide  of  radius  rg.  Both 
the  beam  radius  and  beam  centroid  displacement  are  as¬ 
sumed  to  be  small  in  comparison  to  the  waveguide  radius. 
We  expect  TEn  mode  to  have  the  largest  growth  rate, 
because  its  electric  field  peaks  on  axis. 

The  dispersion  relation  for  t^e  TE-,  t  mode  in  the  pres¬ 
ence  ot  the  beam  is  a  function  of  the  determinant  of  a 

‘FM  Technologies,  Inc.,  Fairfax,  VA 
Manuscript  approved  August  7.  1991 


tri-diagonal  matrix  of  the  form 
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and  K\  =  Kqkqj\f2.  Here,  and  k  are  the  frequency 
and  wavenumber  of  the  TEu  mode,  Kq  -  Qq/v0  is  the 
wavenumber  associated  with  the  quadrupole  field,  where 
=  \e\Bq/-y0m0c,  k\  =  2u/y0ln,  v  ~  IJ\10n,  h 
is  the  electron  beam  current  in  units  of  kilo- Amperes, 
hi  =  A*ii2(Miir7  —  l)*7i  (MuT'g),  and  pnr,  is  the  the  largest 
positive  zero  of  Bessel  function  J[. 

It  can  be  shown  from  Eq.  (2)  that  the  growth  rate 
of  the  instability  is  periodic  in  k.  For  a  given  unstable 
frequency  u>0,  the  unstable  wavenumbers  are  at  all  k  = 
Aso  +  nkq,  where  n  =  0,  ±1, ....  is  an  integer,  and  ktl  is  the 
unstable  wavenumber  associated  with  a  vacuum  waveguide 
mode.  Coupling  to  an  infinite  number  of  modes  may  be 
avoided  in  the  approximation  that  4K%/kq  <C  1-  To  zeroth 
order,  the  approximate  dispersion  relation  is 
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To  first  and  second  order,  we  find 
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respectively. 

II.  Analytical  Results 

Much  insight  can  be  obtained  from  analytical  decom¬ 
position  of  the  approximate  dispersion  relation,  Eq.  (6). 
This  can  give  us  analytical  expressions  for  the  growth  rates 
and  instability  boundaries  in  parameter  space. 

The  dispersion  relation  can  be  rewritten  with  the  cur¬ 
rent  coupling  terms  grouped  together  in  a  term  a, 

w0Mr+1ir_1n  =  ff,  (9) 

where 

Wm{u),k)  =  (tv/c)2  -  {k  +  rnkq)2  -  (10) 

n(u>,fc)  =  -  y  [<4i  +<*li] .  (ii) 

am(w,  k)  =  u/v0  -  (k  +  mfc9),  (12) 

and 

9  =  ifejctfa*  i«li  +  W+iW-x  +  \VnlV_,  +  H'nW'+I] 

-\klK}W0  [o*  jW.,  +  a*  !«•'+,]  .  (13) 

The  polynomial  I]  can  be  rewritten  as 

n(w,  k)  =  {an  -  kq)2A  -  2k,K$a„,  (14) 

where 


Since  is  typically  much  less  than  1,  the  last  term  in 
(14)  can  be  dropped. 

There  are  six  waveguide  modes,  i.e.,  WoW7+iW,_1,  and 
there  are  six  beam  modes.  Four  of  the  beam  modes  in 
(14)  have  the  same  characteristics  as  those  in  the  helical 
quadrupole  case.2  As  in  the  helical  quadrupole  case,  the 
electron  beam  develops  orbit  instability,  when 

(fc2/4-K2)<0,  (16) 

with  or  without  electromagnetic  modes  in  the  waveguide. 
Electromagnetic  instability  exists  at  the  intersection  of  the 

Wr_i(u>,  k)  =  Q  (17a) 

waveguide  me" -  mid  the 

<c  -4-  {-kg/2  -  y^/4  -  A2)  -  ui/vn  -  0  (176) 

beam  mode. 


Following  the  same  procedure  as  outlined  in  Ref.  2, 
analytical  expressions  for  the  spatial  or  temporal  growth 
rates  can  be  obtained.  The  derivation  for  the  spatial 
growth  rate  is  presented  here.  Electromagnetic  instabil¬ 
ity  occurs  at  frequency  ui  which  satisfies 

\J{ujc)2  -  Mn  +kf  =  (w/t>0)  + V2“  \Al/4  ~  Kd-  (18) 
The  dispersion  relation  (9)  can  be  rewritten  as 

{k2  -  k2){(k  +  k„)2  -  k2)((k  -  kq)2  -  k2)(k  -  k2)2 

x[(Jfe  -  k4)2  -  Aifc2][(Jfe  -  k4)2  -  Aik2]  ~  -<r,  (19) 

where  kt)  =  v/(u?/c)2  -  fi2n,  kx  =  w/vn,  k2  =  (u/v..)  -  kq, 
k3  =  (w/vn)  +  kq,  fc4  =  (u)/vo)  +  kq/2,  A kx  =  yjk2/ 4  -  Kj, 

and  A k2  =  y^|/4  +  K\.  Defining  k  =  fe0  4-  kq  +  6k,  the 
imaginary  part  of  6k  is  the  spatial  growth  rate.  The  ana¬ 
lytical  growth  rate  expression  can  be  simplified  to 

6k2  =  a/Da,  (20) 

where  a  =  d\k=k„+k, , 

o  ~  k2{k  -  k3)2(k2  -  k2)((k  +  kg)2  -  k2) 

{{k-krfik-k,)2-*!}  (21) 

and 

Da  ~  -AAk^Ak2  -  A k2)kn(k  -  k2)2 
x(k2  -  *2)((*  +  kq )2  -  (22) 

When  the  more  complete  dispersion  relation  (8)  is 
considered,  growth  rate  at  other  intersections  in  (tv,fc)  of 
waveguide  modes  and  other  beam  modes  will  collapse  to 
an  expression  similar  to  Eq.  (20)  (with  the  wavenumber  k 
shifted  by  ±nkq). 

III.  Examples 

Here,  the  dispersion  relations  (8)  and  (9)  are  solved 
numerically.  We  verified  that  i)  Eq.  (9)  is  a  fair  approxi¬ 
mation  to  (8),  ii)  the  beam  mode  decomposition  of  (9)  is 
accurate,  iii)  the  boundary  of  orbit  instability  is  as  pre¬ 
dicted  by  Eq.  (16)  and  iv)  analytical  growth  rate  expres¬ 
sions  are  in  good  agreement  with  the  results  from  (8). 

Parameters  for  numerical  examples  are  typical  of  a 
high-current,  induction-accelerator  beam,  Ir  =  1  kA  and 
7,  =  7  We  chose  a  quadrupole  gradient  of  Bqkq  =  200 
G/cm,  quadrupole  wavenumber  of  kq  —  0.5  cm-1  and  the 
waveguide  radius  of  rg  =  3  cm.  Figure  1  plots  the  dis¬ 
persion  relation  (9).  The  solid  circle  in  Fig.  1  marks  the 
instability  denoted  by  Eqs.  (17a,b).  Figure  2  has  plots 
of  the  spatial  growth  rate  as  a  function  of  wavenumber, 
k,  for  tv  >0  from  the  improved  dispersion  ration  fS) 
Thv  instability  is  periodic  in  k.  The  growth  rates  of  all 
the  instabilities  are  identical,  except  one.  The  instabilities 
at  ( k  ~  fcn  —  kq  and  k  ~  kn  —  2 kq)  are  from  coupling  to 


T+  2-  The  instabilities  from  coupling  to  T_2  fnot  shown) 
are  at  (k  ~  +  2 kq  and  k  ~  ko  +  3k,,).  The  instability 

at  ( k  ~  k0  -  2 kq),  similar  to  the  instability  marked  by  the 
dashed  circle  in  Fig.  1,  is  erroneous,  and  is  modified  to 
give  the  same  growth  rate  as  that  at  k0  +  kq  as  more  terms 
of  the  full  determinant  (2)  are  included.  The  analytical 
result  from  Eq.  (20)  is  indicated  by  the  x  on  Fig.  2. 


linear  accelerators  because  the  electron  beam  pulse  length 
is  short.  The  instability  group  velocity  is  less  than  the 
speed  of  light  and  propagates  out  of  the  tail  of  the  beam 
pulse  before  it  can  grow  substantially. 

This  work  was  supported  by  the  Defense  Advanced 
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Fig.  2.  Spatial  growth  rates  from  Eq.  (8)  for  w  >  0. 
IV.  Summary  and  Comments 

Electromagnetic  instability  was  found  to  grow  on  an 
intense  beam  in  a  FODO  lattice,  Exact  and  approximate 
dispersion  relations  were  derived.  Analytical  growth  rates 
were  obtained  and  are  in  good  agreement  with  numerical 
solutions  of  the  dispersion  relation.  The  boundary  between 
•'ectromagnetic  and  orbit  instabilities  was  obtained  and 
verified.  The  analysis  will  h»  -vt»nded  to  consider  cases 
which  include  an  additional  nonzero  axial  (solenoidal)  field. 

Note  that  this  instability  has  not  been  observed  in  rf 
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